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^ ■ ABSTRACT 
00 ' 

Q \ Software package for Monte-Carlo simulation of e^e~ exclusive annihilation channels 

' written in the C-|--|- language for Linux/Solaris platforms has been developed. It 

^ ■ incorporates matrix elements for several mechanisms of multipion production in a 

^ . model of consequent two and three-body resonance decays. Possible charge states of 

^ \ intermediate and final particles are accounted automatically under the assumption 

Q \ of isospin conservation. Interference effects can be taken into acccount. Package 

^ ' structure allows adding new matrix elements written in a gauge-invariant form. 
> 

X 



' 1 Introduction 



Simulation of hadron production at low energies {^/s ~ 1 — 3 GeV) is relevant 
for various physical problems. As the cross-section in this region is considered to 
be saturated with intermediate resonances, it is possible to study the properties of 
these resonances through mass, momentum and angular distributions. Simulation 
of different channels allows one to better account for selection rules and interfer- 
ence effects. In particular, the problem of measuring R (total cross-section of e^e~ 
annihilation into hadrons) is a crucial information for contemporary high energy 
physics, especially for the ^^y^ problem jT]. At low energies this value can only 
be obtained from an experiment. Simulation allows to calculate the ratios between 
different channels and estimate the accuracy of measurements. 
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Of other related problems we can highlight the studies of e^e~ Gvr at 
^/s < m-r and (Gtt)^//^ decays allowing to check the vector current conservation 

hypothesis in the Standard Model [21 El, studies of the p" Gtt decay mode as well 
as D- and S-meson decays. 

2 Package features 

At this moment there is no theory, describing the strong interactions reliably at low 
energies. We used a common phenomenological approach, assuming the production 
of a final state in consequent resonance decays: 

• only tree-level diagrams including 2- and 3-body decays, are considered; 

• there may be several interfering mechanisms (e.g. 7* n'^af , p^a) which may 
contribute to different final states (e.g. 2tt~^2tt^ ,n^n^2n^) in accordance with charge 
and strong isospin conservation; 

• possible permutations of final particles are taken into account. 

To simplify adding new matrix elements, we introduced the following requirements: 

• the tensor structures are chosen according to the field transformation properties 
(spin, isospin, parities); 

• expressions for matrix elements are written in gauge-invariant form with the help of 
specialized tensor library (e.g., e^ecTTPp^„^,e^emPa^{)e^u-t&P°'^^e^uapPpi2a{PTT^ - P-k2)p)\ 

• indices' contraction and particles permutations are performed by the software. 
Allowed charge states and relative phases of permutations are based on the strong 
isospin part of the matrix element. Scalar parts of propagators may have arbitrary 
forms, including form-factors and dependence of width on virtuality. In case of the 
vector intial state (7*), the absolute value of the transverse part of hadronic current 
is used as the value of matrix element. Several interfering matrix elements with 
arbitrary complex relative coefficients may be used. 

3 Comparison with experimental data 

By now the package has been used to calculate 4 and 6-pion production cross-sections 
in e"'"e"-collisions. Some distributions of simulated data vs. experimental results jl] 
are given in Fig. ^ (parameters not fitted) . Comparison with analytical calculations 
given in Ref. jH] show a reasonable agreement. 

At the VEPP-2000 collider (the major upgrade of VEPP-2M, BINP, Novosibirsk) Gyr 
production will become possible. Despite the lack of absolute normalization factor 
in matrix elements, the ratio of cross-sections of different final charge states with 
the same production mechanism can be tested for Gtt data. 
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Figure 1: Distribution on dipion invariant mass m^+^o in the channel e'^e aivr 
pTiTT — >• 7r"'"7r~27r° - our simulation (left) and experimental data (right). 

4 Conclusion 

The package for simulation of multipion production has been developed. It includes 
generators for 7* ^ pyr 3%, ai7r,p/o — 47r, ujtttv Stt, and ljtt, pfQ{pp, aa) — > 
6tt production mechanisms. This tool may be useful for many problems involving 
production of multipion final states. 
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